Magnetotransport properties of a polarization-doped three-dimensional electron slab 
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We present evidence of strong Shubnikov-de-Haas magnetoresistance oscillations in a polarization- 
doped degenerate three-dimensional electron slab in an Alj;Gai_a:N semiconductor system. The 
degenerate free carriers are generated by a novel technique by grading a polar alloy semiconductor 
with spatially changing polarization. Analysis of the magnetotransport data enables us to extract an 
effective mass of m* = 0.19mo and a quantum scattering time of Tq = 0.3ps. Analysis of scattering 
processes helps us extract an alloy scattering parameter for the Ala:Gai-a;N material system to be 
Vo = l.SeV. 
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I. INTRODUCTION 

In a recent paperEI, we reported the creation 
of polarization-doped three-dimensional electron slabs 
(3DES) in a graded AlGaN layer. In this work, we 
present magnetotransport studies of the polarization- 
doped 3DES. The 3DES exhibits strong Shubnikov de- 
Haas oscillations. A study of these oscillations reveals 
the effective mass of carriers, the collisional broadening 
due to quantum scattering , and the nature of transport 
in the polarization-doped 3DES. 

In a polar crystal with uniform composition, n = 
V • P = (where P is the polarization vector), and there 
is no bulk polarization charge. The polarization of the 
crystal manifests only as sheet charges at surfaces and in- 
terfaces, owing to the change in polarization across them 
(cr = (Pi — P2) • n, where a is the sheet charge. Pi, P2 
are the polarization vectors on the two sides of the junc- 
tion, and n is the normal to the junction). An alloy 
of spatially constant composition of two polar materials 
with different magnitudes of polarization will also have 
n = V • P = 0. However, if the composition changes spa- 
tially, the divergence of the polarization vector becomes 
non- vanishing. This leads to an immobile polarization 
charge = V • P in a graded alloy region. 

The III-V nitride family of crystals of the wurtzite crys- 
tal structure is among the first material systems that is 
strongly polar (compared to other III-V semiconductors) 
along the c-(OOOl) axis, is semiconducting (most strongly 
polar crystals are insulators and used as dielectrics), its 
binaries (GaN, AIN, InN) have different polarizations, 
and they can be grown expitaxially with sufficiently pure 
crystalline quality. These unique properties of the ma- 
terial system enables us to exploit the concept of po- 
larization bulk doping. The large difference in polar- 



ization between AIN and GaN (both spontaneous and 
piezoelectric)l3 and the good control of epitaxial growth 
of AlGaN alloys prompted its choice for this experiment. 

Figure I shows a schematic of charge control and band 
diagram of the technique of polarization doping that 
we have employed. Also shown in the figure is the 
sample structure we have used. The sample is a Ga- 
face strjticture grown by plasma-induced molecular beam 
epitaxyO on a metal-osganic chemical vapor deposition 
grown semi-insulatingjj GaN on a sapphire, substrate. 
The growth is along the polar c(OOOl) axisH. The top 
lOOnm of the structure is linearly graded AlGaN; the 
composition of Al is changed from 0-30% by controlling 
the aluminum flux by a computer programEJ. This is the 
layer that is polarization doped. p 

The linear grading leads to a fixed chargeCI doping 
= V • P terminated by a opposite sheet charge 
= |P|/e at the surface. The large electric field cre- 
ated by the jellium of fixed polarimtion charge extracts 
electrons from surface donor stateaa to form a mobile 3- 
dimensional electron slab (3DES). The mobile 3DES is 
characterized by a temperature independent carrier con- 
centration for OAK < T < 300K verifying that the car- 
riers are degenerate, and thus distinct from carriers ac- 
tivated from shallow donorsEl. We now proceed to study 
the magnetotransport properties of the 3DES created by 
this novel process. 

II. MAGNETOTRANSPORT 

In the presence of a quantizing magnetic field, the un- 
perturbed 3-dimensional density of states (DOS) go{e) 
undergoes Landau quantization to quasi 1-dimensional 
densitv-pf states and acquires an oscillatory component. 
Dingldl3 showed that inclusion of collisional broaden- 
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FIG. 1: Schematic of charge control showing polarization 
charges and formation of the 3DES. The band diagram shows 
depletion of the 3DES from the surface potential. Also shown 
is the epitaxial layer structure that generates the 3DES. 



ing removes the divergence at the bottom of each 1- 
dimensional subband and damps the DOS oscillation am- 
plitudes exponentially in 1/B. 

As is well known from the theory of the magnetic quan- 
tum effects, this oscillation of the density of states man- 
ifests in oscillations of both the diamagnetic susceptibil- 
ity (manifesting in the de-Haas Van Alphen effect) and 
transport coefficients (manifesting in the Shubnikov de- 
Haas or SdH oscillations). In particular, the transverse 
(B 1. E) rMgnetoresistance Rxx shows oscillations in 
1/B. Kubotll derived the transverse magnetoresistance 
at high magnetic fields using the density-matrix approach 
to solve the transport problem. The expression for R^x 
can be decomposed-into a background part and an oscil- 
latory contribution^ Rxx = Rxx"'' + ^Kf- The back- 
ground term is attributed to sample inhomogeneities and 
disorder. The amplitude {Al of the oscillatory compo- 
nent can be cast in a formllj that is simple to use and 
captures the physical processes reflected in the measured 
magnetoresistance - 



A = 



X 



sinhx 

Dt{T) 



X exp(- 



27rr ^ x(-^)^/^cos(^— 



a 

o 
o 



^ 5 



— I — I — I — I — I — I — I — I — I — I — I — r 

T=400nnK 




2 4 6 8 10 12 
B{J) 



FIG. 2: Magnetotransport measurement data at T = AQOmK. 
The figure also shown as an inset the geometry used for mea- 
suring Rxx, Rxy The Rxx shown is the oscillatory component 
with the background removed. 
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where Ai?"*"^ is the oscillating part of the mag- 
netoresistance with the background removed, ep = 
(nV2TO*)(37r2n3£,£;s)2/3 is the Fermi energy of the 3DES, 



h being the reduced Planck's constant and tildes is the 
carrier density of the 3DES, and ujc — eB/m* is the cy- 
clotron frequency, x = '2,i:'^kBT /hwc is a temperature de- 
pendent dimensionless parameter and F is the coUisional 
broadening energy due to quantum scattering events. 5 
is a phase factor that is unimportant for our study. The 
terms Dt{T) and Dc{B) are the temperature and colli- 
sion damping terms respectively; it is easily seen that in 
the absence of damping of the oscillations due to temper- 
ature (lim^^o Dt{T) = 1) and in the absence of damping 
due to collisions (limr^o -Dc(5) = 1)7 the magnetoresis- 
tance would exhibit a weakly modulated (~ B^^^) co- 
sine oscillations in 1/B. In fact, the two damping terms 
Dt{T), Dc{B) are used as probes to tune the temperature 
and magnetic field independently to extract the effective 
mass and the quantum scattering time. The period of 
the cosine oscillatory term yields the carrier density of 
the 3DES since the period is linked to u^des- Rxy, the 
Hall resistance is linear with B, and should show plateaus 
^j'j at the minima of Rxx when a small number of Landau 
levels are filled. 

For magnetotransport measurements on our 3DES, 
ohmic contacts were formed in a Van-der Pauw geometry 
(Figure 1 inset). The sample was immersed in a ^TJe low- 
temperature cryostat with a base temperature of ZQQmK . 
Magnetic fields in the range OT < B < lAT were applied. 
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Rxx and Rxy was measured as in the geometry depicted 
in the figure using the standard low-frequency lock-in 
technique. 

In Figure II we show a plot of the oscillatory transverse 
magnetoresistance Ai?°^^ and R^y at T — AOOmK plot- 
ted against the applied magnetic field. We have removed 
the background using a FFT filter for plotting Rxx- We 
will briefly describe our FFT process later in this work. 
The Hall mobility determined from the slope of the Rxy 
curve is fin — SOQQcrn^ /V ■ s, which is higher than 77K 
low-field Hall mobility of firrx — 2500cm^ /V ■ s. Also, 
assuming that the 3DES is spread over a thickness d, 
the sheet carrier density of the 3DES is calculated to 
be n^DES X d = I/Rhb = B/eRxy = 7.2 x lO^^/cm^. 
This is consistent with the 77K low-field Hall measured 
sheet density of 7.5 x lO^^/cm^. The spread of the 3DES 
is calculated from a self-consistent Poisson-Schrodinger 
band calculation to be d = 75nm due to of 25nm de- 
pletion of the 3DES from the surface potential. This 
depletion in the graded AlGaN layerJias also been ver- 
ified by capacitance-voltage profiling^. Thus, the Hall 
3-dimensional carrier density is u^des ~ 10^^/cm'^. 



III. ANALYSIS OF MAGNETOTRANSPORT 
DATA 
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FIG. 3: The oscillatory component Ai?°J'^ plotted against 
1/B. The oscillations are periodic with period A(1/_B) = 
0.0294r~^ , and are damped with both increasing temperature 
(different curves), and increasing 1/B. Also shown in the inset 
is a typical FFT power spectrum (at T = 2.5A') showing a 
peak at the fundamental period, and the band-pass window 
used to filter the oscillatory component AR"^"^. 



We begin our study by analyzing the oscillatory com- 
ponent of the transverse magnetoresistance (Figure III). 
For achieving this, we first take the raw Rxx vs B data 
and interpolate it to create an equally spaced N = 2^^ 
size FFT window. We then find the FFT power spec- 
trum. This is repeated for Rxx measured at different tem- 
peratures. A typical FFT power spectrum (at T=2.5K) is 
shown in the inset of Figure HI. There is a clearly resolved 
peak at the fundamental oscillation period Bq = 34.0ir. 
A band pass filter [fpass =28-150T] is then employed to 
remove the background component. The resulting AR°^^ 
for various temperatures OAK < T < 9.5K is plotted 
against 1/B in Figure HI. As is clear from the plot, the 
period of oscillations is A(-i) = 0.0294T-1 = I/Bq. The 
oscillations are strongly damped with increasing 1/B as 
well as with increasing temperature, as predicted by the 
theory (Equation 1). 



A. Carrier concentration 

First, we observe from Equation 1 that the period 
A(-g) is linked to the carrier density of the 3DEG by 
the relation A{l/B) = eh/m*SF = ^{Sn'^nsDEs)''^^^ ■ 
From the plot, the period A{l/B) = 0.0294T-^ yields a 
direct measurement of the 3-dimensional carrier concen- 
tration n^^^ — 1.1 X lO^^cm"'^. Thus, the carrier den- 
sity measured from the quantum oscillations is close to 
the carrier density measured by classical Hall technique 



B. Effective mass 

The second important property of the 3DES that is 
accessible from the SdH oscllations is the effective mass 
of the electrons. This can be calculated from a controlled 
temperature damping of the oscillation amplitudes, i.e., 
by tuning the term Dt(T) (which is the opljy tempera- 
ture dependent term). The approximationllj sinh(x) ~ 
exp{x)/2 converts Equation 1 into 



ln(^)^C-(^^!^)T 
^t' ^ ehB ' 



(2) 



where the constant C absorbs all terms independent of 
temperature. Thus, when ln(y) is plotted against T for 
a fixed magnetic field, it results in a straight line with 



slope S = 



ehB 



Since the effective mass is the only 
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unknown, it can be calculated. We calculate the effective 
mass by employing this procedure. 

We choose the amplitude maxima at three different 
magnetic fields B = 8.9T, 10. ST, 12.0r for extracting ef- 
fective mass of the 3DES. These maxima are indicated 
by arrows in Figure 3. The plots of hi(A/T) against 
temperature {2.5K < T < 14.51^) for these magnetic 
fields are shown in Figure IV(a). The slopes of the 
plots at B = 8.9r, 10. 5r, 12.0r yield effective masses 
TO* = 0.189, 0.197, 0.189too respectively. The accepted 
effective mass of electrons in bulk GaN is to* = 0.2too. 
SdH measurement of electron effective mass in bulk GaN 
is difficult since the activated carriers- in donor-doped 
samples freeze out at low temperaturestJ. Since the elec- 
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trons in our 3DES are degenerate, they resist freezeout 
effects and remain independent of temperature. How- 
ever, the 3DES electron is not in bulk GaN but in an 
alloy with a changing composition; we expect the effec- 
tive mass to be a spatial average. The effective mass 
of electrons in AIN is predictedEJ to be m* = 0.32mo; 
a linear interpolation gives an expected effective mass 
of TO* = 0.21 Too at an average alloy composition of 
(x) = (0.75 X 0.3)/2 = 0.11 for the mobile 3DES. 
This is slightly higher than our measured average value 
to*^ = 0.19too. The result is reasonable within hmits of 
experimental error. We also note that the effective mass 
of electrons in two-dimensional electron gases (2DEGs) 
at AlGaN/GaN heterojunctiati&jbas been measured by 
SdH oscillation techniqueli3,EjjiIl,t3 and reported to be 
various values around m* ~ 0.2too, which is close to our 
measured value. 



C. Scattering times 
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FIG. 4: Effective mass plot and Dingle plot of the 3DES. For 
effective mass, ln{A/T) is plotted for three values of magnetic 
fields {B = 8.9r, 10.5T, 12.0T) against temperature. The 
slope yields the effective masses 0.189mo, 0.197mo, 0.189mo 
respectively. The Dingle plot yields a quantum scattering 
time of Tq = 0.27ps at T = 400mK and an average value of 
T,"" = 0.3ps. 



The third important parameter of the 3DES that is 
measured from the SdH oscillations is the coUisional 
broadening energy T (due to Dingle). This term is a 
measure of the smearing of the delta-function discontinu- 
ities in the DOS due to quantum scattering events, and 
it appears as the imaginary part of the single-particle self 
energy function. CoUisional broadening energy is linked 
to the quantum scattering time Tq and the Dingle tem- 
perature Td by the relation T = h/2Tq = "nkBTo- This 
quantity is experimentally accessible from a controlled 
Landau damping of the oscillation amplitudes with \/B 
at a fixed temperature; in other words, by tuning Dc{B). 
Equation 1 can be cast in the form 



ksT — 0.09mey at T = lA', thus satisfying the condi- 
tions required for clear Shubnikov de-Haas oscillations. 

Whilst Tq (the quantum lifetime of an electron in a 
magnetic quantum state) is determined by all scattering 
events, the transport (classical, or momentum) lifetime 
Tc is weighted by a scattering angle factor (1 — cos 9), 
where 9 is the angular deviation in the scattering event. 
Thus, 



— = / P{9)d9 

To J 



(4) 
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A* 



(i^)*A(r) 



eTQ 



B 



(3) 



for extracting Tq and the related quantities F, Tq. Here 
A* are the extrema points of the damped oscillations, 
forming the exponentially decaying envelope. Equation 
3 suggests that a plot of the natural log of the left side 
quantity against 1/B ('Dingle plot')should result in a 
straight Hue whose slope is —■Km*/eTq. Since we have 
already measured the effective mass to*, we can extract 
Tq from the slope. 

Figure IV(b) shows the Dingle plot for T = 400m AT, 
yielding a = 0.27ps. An averaging of the quantum 
scattering times over a range of low temperatures yields 
a value r^" = 0.3ps. The quantum scattering time does 
not show any discernible trend with temperature in this 
range. We calculate the corresponding level broaden- 
ing F = l.lmeV and Dingle temperature To — AK. 
We mention in passing that the Landau level separation 
at -B = lOT is hioc = ^-imeV , sufficiently larger than 
both the measured coUisional broadening of the Lan- 
dau levels (F = lAmeV) and the thermal broadening 



as opposed to the momentum scattering time, which 
is averaged by the angular contribution for small angle 
scattering 



— = j - cos 



(5) 



here P{9) is the transition probability between the 
initial and the scattered states determined by Fermi's 
golden rule, ft can be seen that for isotropic scatter- 
ing events with no angular preference {P{9) is indepen- 
dent of 6), the quantum and classical scattering times 
are the same Tc/t^ = 1. If the dominant scattering pro- 
cess has a strong angle dependence, the ratio is much 
larger than unity. This fact has been utilized to identify 
the dominant scattering mechanism in modulation-doped 
AlGaAs/GaAs two-dimensional electron gasest^l. 

The low-temperature Hall mobility gives us a direct 
measurement of the classical scattering time for the 3DES 
via the relation = eTc/m*; the value is t^ = 0.34ps. 
Within limits of experimental error, the ratio Tc/t^ ~ 1, 
i.e., is close to unity. 



5 



D. Scattering mechanisms 

As Hsu and Walukiewicz have showr@, one has to ex- 
ercise caution before declaring that the dominant scatter- 
ing process in the 3DES is of short-range nature. They 
show that it is possible to have the ratio close to unity 
if there are different scattering mechanisms responsible 
for the quantum and classical lifetimes. At the low 
temperatures at which we measure the scattering time, 
both acoustic and optical phonon modes of scattering 
are frozen out and the charges from any unintentional 
background donors arc frozen onto the dopor sites, ren- 
dering them neutral. Size effect scatterin^^ that occurs 
if the width of the 3DES is much less than the mean- 
free path of electrons is negligible since our 3DES has a 
mean free path A = hkpfj,/e ~ 60nm whereas the width 
of the 3DES is do « 75nm. The chief scattering mecha- 
nisms that can affect mobility are alloy disorder scatter- 
ing (since the 3DES is in a graded alloy), charged disloca- 
tion scattering (owing to the high density of dislocations 
Ndisi ~ lO^cm"^), and ionized impurity scattering from 
the remote donors at the-surface states. 

Hsu and WalukiewicaH3 show that remote ionized im- 
purity scattering strongly favors small angle scattering, 
thus causing the ratio Tc/t^ ^ 1. Since Tc/t^ « 1 for 
our 3DES, remote ionized impurity scattering is unim- 
portant. 

The ratio of classical to quantum scattering times due 
to chaijged dislocation scattering has not been found yet. 
LookE3 recently applied the classical scattering time due 
to charged dislocations and applied it to explain the mo- 
bility of bulk doped GaN samples. The ratio classical 
to quantum scattering times due to charged^dislocation 
scattering in a degenerate 3DES is given byCJ 

4^ = 1 + 2fc2 (6) 

where kp — {Stt^ti^ijesY^^ is the Fermi wavevector 
and A^^ — 2t£p /ie^ tildes is the Thomas-Fermi screen- 
ing length for a degenerate 3DES. The ratio, which de- 
pends only on the 3D carrier density for a degenerate 
3DES, is evidently greater than unity; for our 3DES, 
the ratio is 2.3, which is larger than what we observe. 
So dislocation scattering is not the dominant scattering 
mechanism for determining the quantum scattering rate. 
The dislocation scattering dominated classical scatter- 
ing rates is strongly dependent on the amount of charge 
in the dislocation cores, which-4p low for our 3DES due 
to the lack of shallow dopantsa. Thus, we exclude dis- 
location scattering to be the most important scattering 
mechanism for classical scattering rate. 

So we converge on alloy scattering as the dominant 
scattering mechanism for low temperatures. The quan- 
tum and classical scattering rates for alloy disorder scat- 
tering are the same since the scattering potential Vo is of 
a short range nature, which makes the scattering process 
angle independent. 



The scattering rate due to alloy ^disorder with a short 
range potential Vb has been showr£3 to be given by 

1 2lT 

— = —V^n{x)x{l - x)gMek) (7) 

Talloyik) n 

where i^o{x) is alloy composition-dependent volume of 
the unit cell over which the alloy scattering potential Vq 
is effective, and x is the alloy composition. <73d(£) is 
the 3-dimensional DOS. Since the 3DES is highly de- 
generate, only carriers with wavevectors very close to 
k = kp contribute to transport at low temperatures; 
hence, the scattering rate is evaluated at the Fermi en- 
e^gy '^ailoy = '^aiioyi^p)- Bcsides, the alloy is graded, and 
Matheissen's rule which is accurate for low temperature 
transport analysis is used for a spatial averaging of the 
scattering rate 

Kiloy) ^ ^ f ° r-l^ix)dx (8) 

where xq — 0.225 is the alloy composition experienced 
by 3DES electrons at the edge of the depletion region. 
Using this simple result, we conclude that to achieve a 
low-temperature transport mobility of SOOOcnr^ /V ■ s, an 
alloy scattering potential of Vq = l.SeV is necessary. 

It is important to point out that the alloy scattering 
potential as defined in the model for alloy disorder is at 
best a fitting parameter. Due to the lack of experimen- 
tal values, it is common practice to assume the scatter- 
ing potential to be the conduction band offset between 
the binaries forming the alloy (Vb — AEc — 2.1eV for 
AlN,GaN). We note that with an alloy scattering po- 
tential of Vb = 2.1eV, the calculated mobility is lower 
(« 2000cm^ /V ■ s) than the measured value. Besides, 
our 3DES mobility is dominated by alloy scattering and 
all other scattering mechanisms are removed by measur- 
ing the mobility at low temperatures, making it a clean 
measurement of the alloy scattering potential. Our re- 
port presents the first direct measurement of the alloy 
scattering potential in the AlGaN material system. 

IV. SUMMARY 

We finally summarize our findings in this work (see Ta- 
ble I). By exploiting the polarization charges in the Al- 
GaN/GaN semiconductor system, we are able to create 
a 3DES without intentional doping. The mobile carriers 
in the 3DES are degenerate and exhibit a high mobil- 
ity even to the lowest temperatures. The lack of carrier 
freezeout enables us to observe Shubnikov de-Haas oscil- 
lations in magnetotransport measurements of the 3DES. 
The oscillations reveal several important facts about the 
3DES. First, the temperature damping of oscillations re- 
veals the effective mass of electrons to be very close to 
that in bulk GaN {m* = 0.19mo). Next, the quantum 
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scattering time of electrons in the 3DES is found from 
the Dingle plot to be = 0.3ps. The ratio of the clas- 
sical (momentum) scattering time to the classical (mo- 
mentum) scattering time is found to be close to unity 
Tc/tq ~ 1. The ratio suggests predominantly short-range 
scattering dominating transport properties at low tem- 
peratures. This scattering mechanism is identified to be 
alloy scattering. This lets us extract another valuable 
parameter, the alloy scattering potential in Al^rGai-^^N 
to be Vo = l.SeV. 
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TABLE I: Summary of constants used and results extracted from magnetoresistance measurements 



Quantity 


Symbol 


Magnitude 


Unit 


Relative dielectric constant 




8.9(GaN), 8.5(A1N) 




Lattice constant 


ao 


o ^ /A / /^i '\T\ o -t -t c\ / Al"\T\ 

3.189(GaN), 3.112(A1N) 


A 


Lattice constant 


Co 


5.185(GaN), 4.982(A1N) 


A 


Effective mass 


m 


n 1 n 

u.iy 


mo 


Quantum scattering time 


Tq 


0.3 


ps 


Transport scattering time 


Tm 


0.34 


ps 


3DES density 




1.1 X 10^* 


cm~^ 


3DES Hall mobility (IK) 




3000 


err? jV ■ s 


Alloy scattering potential 


Vo 


1.8 


eV 



